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Abstract — In this paper, we consider the distributed robust 
filtering problem, where estimator design is based on a set 
of coupled linear matrix inequalities (LMIs). We separate the 
problem and show that the method of multipliers can be applied 
to obtain a solution efficiently and in a decentralized fashion, 
i.e. all local estimators can compute their filter gains locally, 
with communications restricted to their neighbors. 

I. Introduction 

Estimator design has been an essential part of controller 
design ever since the development of state-space based 
controllers. A milestone was laid by the Kalman Filter in 
1960 [1], 

While in the classical estimator design one estimator is 
used for the entire system, distributed estimators have gained 
attention since a distributed Kalman Filter was presented in 
[2], [3]. In a distributed estimator setup, multiple estimators 
create an estimate of the system’s state, either individually 
[4] or cooperatively. In the latter case, even when every single 
estimator may be able to obtain an estimate of the state 
on its own, cooperation reduces the effects of model and 
measurement disturbances [5]. Also, the situations are not 
uncommon where individual estimators are unable to obtain 
an estimate of the state on their own and cooperation be¬ 
comes an essential prerequisite [6], [7]. The node estimators 
may even not have a model of the full system, but only know 
a part of the system [8]. 

However, even though the setup consists of distributed 
estimation units without a central coordinator, in many 
known approaches the design process itself requires a central 
coordination unit. In some practical application examples, 
where the design process can be done offline, this may 
not be a significant drawback. On the other hand, in many 
applications especially those involving disttibuted sensor net¬ 
works with varying communication topology, a centtalized 
computation of observer parameters represents a severe lim¬ 
itation. Practicality of a distributed system demands that the 
estimator design process is to be carried out in a distributed 

This work was supported by the German Research Foundation (DFG) 
through the Cluster of Excellence in Simulation Technology (EXC 310/1) 
at the University of Stuttgart and the Australian Research Council under 
Discovery Projects funding scheme (Project DP120102152). 

J. Wu and F. Allgower are with the Institute for Systems Theory and 
Automatic Control, University of Stuttgart, Stuttgart, Germany. Email: 
{jingbo . wu, allgower}@ ist. uni-stuttgart. de. L. Li is with 
School of Electrical, Mechanical and Mechatronic Systems, University of 
Technology, Sydney, Australia. Email: Li.Li@uts.edu.au; V. Ugri¬ 
novskii is with the School of Information Technology and Electrical En¬ 
gineering, University of New South Wales at the Australian Defence Force 
Academy, Canberra, Australia, Email: v. ougrinovski@adf a . edu . au. 
Part of this work was done during L. Li’s visit to the School of Information 
Technology and Electrical Engineering, University of New South Wales at 
the Australian Defence Force Academy. 


manner as well. If network needs to adapt to some changes, 
such as a change in the plant or change in the network 
structure, this allows each node to reconfigure using only 
local communications and computation only. 

In this paper, we provide a complete analysis of one dis¬ 
tributed estimation problem where such a distributed design 
scheme is possible. Specifically, we adopt the setup from [6] 
concerned with the problem of distributed estimation with 
Mio consensus performance. As a matter of fact, in [6] a 
gradient-descent-type algorithm was proposed that can be 
used to calculate the filter gains in a distributed manner. 
Although the proposed gradient type algorithm demonstrated 
a possibility of computing the estimator parameters in prin¬ 
ciple, a practical application of that algorithm is hindered 
due to slow convergence observed even in low dimensional 
examples. Also, implementation of the decentralized design 
scheme proposed in [6] requires bidirectional communica¬ 
tions between the network nodes, which essentially requires 
the communication graph to be undirected for the purpose 
of the estimator design. In this paper, we address the prob¬ 
lem of designing distributed estimators by using distributed 
optimization methods presented in [9],[10]. Distributed opti¬ 
mization methods are widely applied in networked systems, 
see e.g. [11],[12],[13]. The contribution of this paper is to 
show that the problem of designing distributed estimators is 
amenable to the methodology of distributed optimization as 
well. Although the design scheme is proposed for the specific 
class of algorithms in [6], it illustrates all the steps neces¬ 
sary to devise similar design schemes for other distributed 
estimation algorithms and distributed optimization subject to 
LMI-constraints in general. 

The rest of the paper is organized as follows: We first 
introduce the notation and some preliminaries on graph 
theory. Then, we revisit some essential results published in 
[6] and discuss there limitations with respect to numerical 
optimization. Section III is dedicated for introducing the 
proposed optimization scheme. In Section IV, we give a 
mathematical example, and Section V concludes the paper. 

II. Preliminaries and Background 

In this section, we introduce the basic definitions and 
results which our main results will build on. 

A. Notation 

Let P be a symmetric matrix. If P is positive definite, it is 
denoted P >0, and we write P < 0, if P is negative definite. 
0 denotes a matrix of suitable dimension, with all entries 
equal 0. Moreover, for vectors x G we use the Euclidean 
vector norm ||x|| = sjxJx and the weighted vector norm 


||;ii:||/> = \/Px for symmetric matrices P > 0. For matrices 
A G we use the Frobenius norm ||A|| = ^ytr{A~^A) and 

the induced norm ||A ||2 = suP||^||^q(||Ajc||/||;ii:||). ^2 denotes 
the Lebesgue space of R"—valued vector-functions z(-), 
defined on the time interval [ 0 ,oo) with the norm ||z ||2 = 
V /(T The vectorization vec(-) maps any matrix 

A G jQ jjjg n • m-dimensional vector vec(A) formed by 

the stacked columns of A. 

B. Communication graphs 

In this section we summarize some notation from the graph 
theory. We use directed, unweighted graphs ^ to 

describe the communication topology between the individual 
agents. 'P = {vi,...,VAr} is the set of vertices, where vt&y 
represents the k-th agent. S' CY x "P are the sets of edges, 
which model the information flow, i.e. the k-th agent receives 
information from agent j if and only if {vj,Vk) G S. The set 
of vertices that agent k receives information from is called 
the neighborhood of agent k, which is denoted by ^ = {j ■ 
{vj,Vk) G S}. The set of vertices that receive information 
from agent k is called the out-neighborhood of agent k, which 
is denoted by = {j '■ {vkYi) ^ Sj. The in-degree and 
out-degree qi^ of a vertex k is defined as the number of edges 
in S, which have Vk as their head and tail, respectively. 

C. Distributed ML estimation and LMI conditions 

We now present the underlying distributed estimation 
problem from [ 6 ]. It involves LMI design conditions, which 
are the main object of interest in the paper. Our main 
objective is to show that these LMI design conditions are 
amendable to a distributed solution by using the multiplier 
method (c.f. [9]) . 

The distributed estimation problem with M,o consensus of 
estimates posed in [ 6 ] involves estimation of the state of the 
uncertain LTI system described by the differential equation 

x = Ax + B^[t), ( 1 ) 

where x G K" is the system state variable to be estimated and 
^ (f) G M 2 is a disturbance function. N estimators are applied, 
each estimator receives a -dimensional measurement 

yk = CkX-\-Di^^-\-Diir\if^{t). ( 2 ) 

In (EJ, Tjj.(f)G M 2 represents the measurement uncertainty of 
the local estimator k. In particular, it is assumed that E^. = 
DkDj ADkO] > 0 . 

Remark 1: The assumption that > 0 is a standard tech¬ 
nical assumption made in nonsingular M„ control problems 
[14]. It is obviously satisfied in the case when all mea¬ 
surements are affected by disturbances, which is evidently 
satisfied in practical applications. This assumption is later 
used to guarantee boundedness of the solution set. 

The estimators form a network of interconnected Moo 
filters of the form 

Xk =Axk + Lk{yk - Ckh) +KkY^ {xj - 4) n-. 


with initial condition 4(0) = 0. Here the matrices 4 G 
and Kk G are the filter gains to be designed. 

As it can be seen in Q, the estimators are distributed, i.e. 
the local estimators create an estimation of the system’s state 
X, solely based on the local output yk and communication 
with neighbouring estimators. The problem in [ 6 ] was to 
determine estimator gains Lk, Kk in (l3]l to satisfy natural 
internal stability and Moo gain conditions. To introduce these 
conditions, define the local estimator error as =x — Xk, and 
the estimator disagreement function is defined as 

= L (4) 

k=ljeJi 

where x= [x]",and e = [ej^...^eJ^Y■ The estimator 
design problem is concerned with achieving the following 
properties; 

(i) In the absence of model and measurement disturbances 
(i.e., when rjk = 0 ), the estimation errors decay so 
that ek ^0 asymptotically for all k= 

(ii) The estimators 0 provide guaranteed Mo performance 

JY^{x{t))dt 

1 TV / , TV \ (5) 

for some positive definite matrix P, some 7 > 0, and 
performance index 7 > 0 . 

Property (ii) requires both the local estimation errors 
and the estimator disagreement to be bounded with respect 
to the disturbances in an sense. As shown in [ 6 ], 
LMI-conditions can be found, where the solution delivers 
estimator gains sufficient for solving the above problem. To 
present these LMI conditions, define the matrices 

Ak =A + Ukl — BDJEj^ 'Cjt, 

Qk = 4 rAi +AjXk — cjEj^ ^Ck + P {pk qk)I, 
h=[B{I-DjEYDk) -BDJeY'^, 

where Xk G R"^" is a symmetric, positive definite matrix and 
dk^P are positive parameters. For the remainder of this paper, 
we will make two assumptions on the system class. 

Assumption 1 The communication graph is connected 
and balanced, i.e. qk = pt for all k = 1 , ...,A. 

Assumption 2 For all k= 1,...,A, the tuple {Ak,Bk) is 
controllable. 

The LMIs used for designing the estimator gains are 
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where Xk,Fk, and j3 are the variables, p > 0 is a constant 
parameter, and ^ = {jii ■•■1 Jpj,}- formulate a 

variation of the main result from [ 6 ]. 


Proposition 1: Suppose the interconnection graph and 
the parameters ak> Q, k=\,...^N, are such that the set 

r = {j3 > 0 ; (| 6 ]l-(l 8 ]l are feasible for k = 1 , (9) 


is not empty. For any j3 G F, one solution to the distributed 
estimation problem under consideration, with 7 = is given 
by the network of estimators Q in which 

Kk=X^^Fk and Lk = {X^^Cj + BDI)E^\ ( 10 ) 

where Xk and Fk,k = 1 , belong to the feasibility set of 
dU - (l 8 ]l, corresponding to this particular value of j3. The 
weighting matrix P in Q is given by P = (1 /N)Ylk=\Xk- 

Remark 2: Assumption 1 is a restriction toward the class 
of communication graphs, which is made in order to ensure 
that the well-known average consensus algorithm is appli¬ 
cable. Assumption 2 is used to ensure boundedness of the 
feasible sets. It is not restrictive, as it represents the worst 
case of disturbance, and if not satisfied, small hypothetical 
disturbances can be added to the system description, i.e. 
additional columns to B.Dk and Dk- Furthermore, note that 
the tuple {Ak,Ck) is not required to be detectable. 

Since the LMIs (l 6 ]l, Q, dDl are coupled, they may be 
solved in a centralized manner as the optimization problem 

min (-j3) 

( 11 ) 

subject to ([6l),(|7]i,®,k = 1 ,...,A, 

where the resulting matrices Xk,Fk deliver the estimator gains 
Lk,Kk according to ( fTOl l. In the next section we will explore 
the separation of the problem and parallel computation in 
order to solve the problem in a distributed manner. 


III. Distributed Calculation of Filter Gains 

Parallel and distributed computation is thoroughly dis¬ 
cussed e.g. in [9], and in this section, we use some of 
the methods presented in Section 3 in [9] to calculate 
our estimator gains in a distributed fashion. Solving the 
optimization problem (fTTI) can be formulated as a separable 
problem by defining local representations of the solution 


variables, Xj, and for all k = 1,... ,N and j = k,j\, 

The tuple of local variables is denoted by 

7l=(Fl,/3^X,^4,...,4 ), (12) 

where the upper index k denotes the representation of a 
variable used by estimator k and all Xj are symmetric, 
positive definite matrices, and > 0 . 


Problem 1: Find an iterative algorithm, which creates a 
sequence Yk{t),t G N, such that local representations of the 
variables converge in the sense that 

Iim(j3''(0-i3'^(0) =0, (13) 

for all ki,k 2 = 1,...,A^ and 

lim(x;‘(f)-^'^(0)=0, (14) 

for all j = and ki,k 2 G j. All iterations Yk{t) 

shall satisfy the LMIs (| 6 ll-(l 8 ]l when setting j3 = P'‘,Xk = 

xKX-k =XK, ...,X± =XK . Furthermore, the iteration steps 
h A Jpk Api, ^ 

of the local variables Yk[t -f 1 ) shall be calculated in a 
distributed fashion, i.e. interaction with the neighbors j G 
only. 

As a first step, in order to ensure that both (HB and 
Problem 1 are well-posed, we establish a statement about 
the boundedness of the feasible set of the LMIs (| 6 ]l-®. The 
proof of this theorem will later be used in order to ensure 
that solutions of local optimizations are always attainable. 

Theorem 1: Suppose the pairs {Ak,Bk) are controllable. 
Then, for any p > 0, the feasible set 

Tl = {{p,Xk,Fk,k=l,...,N)\ 

dH),©,® hold for k= 1 ,...,A} 


is bounded. 

Proof: Suppose {l5,Xk,Fk,k = G Tl. Using 

the Schur complement, it follows from ®, ® that for an 
arbitrary > 0 , 

XkAk +AlXk - CjE^^Ck + I5{pk + qk)I 

—PkFk — PkFk +XkBkBlXk 

+TkiFiJXf^Fk-pXk) < 0 . (16) 


Completing the squares on the left-hand side yields 

XkAk +AlXk - CjEf^Ck + ^{pk + qk)I 
+rk{Fk - —XkVXf\Fk - ^Xk) 

-Tki^+p)Xk + XkBkBlXk <0. (17) 

^k 


Hence, we conclude that {l5,Xk) satisfy the Riccati inequality 


Xk{Ak- 
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2Xk 
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2Xk 


I?Xk 


-ClE,^^Ck + ^{pk + qk)I + XkBkBlXk < 0 . (18) 











After pre- and post-multiplying ( fTSl l by X, ^ ( fTSl ) reduces to 
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,2 I - 2 , 


,7 Pk + ^kP , ^-\,7 + 

-- 

< 0(19) 


Associated with this Riccati inequality, consider the Ric- 
cati equation 

(Ak - + Z,(A, - Pk + hP jy 


2Xk 


2Xk 


1 


-Zk{c; E,;'Ck- -ipk + qk)I)Zk+BkBl =0 (20) 


and dehne 


7 ° = inf 


7 > 0 : equation (l20l i has a nonnegative- 
dehnite solution 


( 21 ) 


From the JZL control theory [14, Theorems 4.8 and 9.7], 
it is known that the set whose inhmum determines 7 ° is 

nonempty if the pair {A^— '‘ 2 A is detectable and 

^22 * 

the pair (A^. — is stabilizable. Note that by the 

condition of the theorem, the pair {Ak,Bk) is controllable; this 
implies the stabilizability of (A^— Now, let us 

choose Ti; > 0 such that all unstable unobservable modes of 
the matrix pair (Ak,Ck) lie in the region Re i < . This 

will guarantee that the pair (A — is detectable. 

Thus, we conclude that 7 ° < 00 . 

The feasibility of the Riccati inequality JT^ also implies 
that the following state-feedback control problem in¬ 
volving the system 

2 2 

x={Ak -— - — I)^x + Cju-\-{pk + qk)^^'^w^ (22) 

2Tk 
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and the performance criterion 


f \\zkfdt<^[ \\wfdt VweL 2 , (.x(O) =0),(23) 
Jo p JO 

has a solution. Indeed, it follows from ([T9] l that 
{Ak - P' + ^P l-X,'cjE-^Ck)X-' 


2Xk 


+Xj^ (Ak - 


2Xk 


k '-k ^k '-kj^k 

I-X^^CjE^^CkV 


127 Pk'^'^kP T v-lr’Tc’-lr’ ^T 


+X^ 1 {CjE^^Ck + ^{pk + qk)I)Xk ‘ + BkB'k < 0. (24) 


Since ‘ > 0 and 


is a strict inequality, the matrix 




(25) 


is Hurwitz. Thus, the closed loop system consisting of the 
system (l 22 l l with w = 0 and the state-feedback controller 

u = -E^'^CkXj^'^x 


is exponentially stable. Also, using the completion of 
squares, it is easy to show from ( [T^ that the above controller 


guarantees the attenuation property (|2^ . Since the 
pair {Ak,Bk) is controllable, these observations guarantee 
that the Riccati equation (l20l i with 7 = p has a unique 
nonnegative dehnite stabilizing solution Zk (e.g., see [15, 
Theorem 3.2.2]). Thus, i > 7°- Furthermore, since (Ak,Bk) 
is assumed to be controllable, Zk>0 and is invertible. 

From Theorem 4.8 in [14], we know that 7 ° > 0. These 
observations imply that j3 < 7 °^^ Also, using the relation¬ 
ship between solutions to the Riccati equation (i20l l and the 
corresponding Riccati inequality ([T9] l [16, Lemma 8.1], it 
follows that Xk <Z^^. 

This discussion leads us to conclude that there exist upper 
bounds on feasible j3 and ||Z^||. Indeed, 7 ° andZ^ are dehned 
using the conditions involving the properties of the matrices 
Ak, Ck and Bk and the constants p, pk- Hence, these constant 
and the matrix are not dependent on the choice of the feasible 
j3 and Xk- 

It remains to show that there is an upper bound on the 
feasible Fk as well. Using the Schur complement, ill is 
equivalent to F^Xj^^Fk < pXk- This further implies 

X-'^^FjlXi;^FkX-^/^<pI 

tr(x,^/^F,lX^^FkX-^/^)<np 

< v^. 

For the Frobenius-norm of Fk, we can now conclude 


which is bounded due to boundedness of Xk- ■ 

Now, the decoupled version of the LMI conditions (| 6 ]l, (l7]i, 
ill is proposed as 
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(28) 


with Ql = X^Ak + Ajxl - CjE-^Ck + liHPk + qk)I- Note 
that the LMI-conditions are formulated as non-strict inequal¬ 
ities, but with additional parameter 5 > 0. However, as 5 can 
be chosen arbitrarily small, it introduces no conservativeness. 

We denote the feasible set of the k-th group of the LMIs 
as Elk = {Yk I (l26l l. (l27l) . (i28T l hold true}. Then, the separable 





























convex program can be written as 

N 

minimize(— ^ j3^) 

tr=l 

subject to Yk G — P 

X^=X,, X],=X., 

J\ J\ 

for every k= Here, ji^Xj,] — are additional 

variables that are needed to make the problem separable. 

Remark 3: The optimization problem ( l29l l can be varied 
in the way that for a given performance parameter j3 > 0, 
filter gains for © are to be found. Then, ( |29] | turns to 
a pure feasibility problem witlmut optimization objective, 
and therefore, the variables and their iterations in the 

following algorithm can be omitted. 


(29) 


^x% =x. 


jpi, 


jpk 


The dual problem has the form 

maximize q(A', (30) 


where = (A^,A^,A\,...,A^^ ) for k = 1,...,A is the 

h hk 

suitable tuple of Lagrange multipliers and the dual function 
q(-) is defined as 


q(A',...,A'^)= inf (31) 

YkeQ.k,k=l,...,N 


L(-) is the augmented Lagrangian function (cf. [9]) 


L(T,,A'^) = f (-^^ + - j3^) + ^ 1^ - 

k=l^ / z 

+ E E Or(AY{Xj-x';)) + ^\\Xj-x^\\A 

k=i jeJiuk^ ^ ^ ^ 

(32) 

with design parameter c > 0. The optimization problem (l30l l 
can now be solved iteratively with Algorithm [T] which is 
initialized with kjt(O) G fljt, A > 0 and symmetric A^>0. 


Remark 4: Out of the three steps in Algorithm 1, clearly 
2) and 3) can be run in parallel by the individual estimators 
separately. Calculation of Step 1 of Algorithm 1 requires 
the evaluation of the mean value, which can be done in a 
distributed manner by applying a consensus algorithm. Under 
Assumption 1, average consensus algorithms can be used to 
calculate j3(f + 1). In particular, discrete time algorithms are 
preferable to keep the concept of an iterative algorithm [17] 
and algorithm which converge in finite-time are useful to 
ensure exact convergence [18], [19]. 

For the calculation of Xk{t + \) in the case of undirected 
graphs, only two steps are needed; All neighbors j G 
pass their Xj^{t) and A^(f) to estimator k. Then, estimator 
k calculates Xk{t + l) and passes it back to its neighbors. 
The calculation of Xk{t + 1) in the case of directed graphs is 
more demanding with respect to the graph topology; Usual 
average consensus algorithms can be applied when for every 
k— 1,..., A, the subgraph 'i^k induced by node k and its out- 
neighborhood ^k^ is a balanced graph. This however can be 
relaxed by adding additional variables XjJ ^ to Yk and 


Algorithm 1: Calculation of iteration step t-rl 

1) Set the fusion variables for j = 1,...,A 

iV iVC 

U'+i) = 2 E E A‘«) 

‘iJ ke^j ke^j 

2) Calculate the new variables for k = 1,...,A 

Yk{t + 1) 

=ar^mm + 1 )_^ 


i:|2 




l)-X^f 


jeJiiJk 

3) Set the Lagrange variables for k = 1,...,A 

A^(f + 1) = A^(f) + c(^(f + 1) - + 1)) 

and for all k = G jVk 


A'j{t + 1) — A‘j{t) +c{Xj{t -f l)—Xj(t + 1)) 


adding Xj = Xj as equality constraint. For instance, if for 
all k^l,...,N Yk = iFk,li'‘,X^,...,X^), then Xkit + l),k = 
1, ...,A can be calculated under Assumption 1 using average 
consensus. This will later be demonstrated in the numerical 
example. 

In order to show the convergence of Algorithm 1, two 
lemmas need to be introduced. 


Lemma 1: The Lagrangian (i32l l can be written in terms 
of the vectorized variables, i.e. 

L(T,,A'^) = £ (-j3^ + A'^(^- + £|^- 

/t=i ^ / z 


+ E E (veciA)yvec{Xj-X^) + U\veciXj-X^)f) 

k=lje^kk>k ^ 


Proof: We have the equalities 


(33) 


tr{A^B) = '^Y^AjiBji = vec{Ay vec{B) 

‘ j 

llAf = rr(A^A)=i:^A^,= ||v.c(A)f. 
i j 


This Lemma shows, that we can recast the problem into a 
problem of a standard form defined on a finite dimensional 
vector space. 

Lemma 2: For fixed ...,Ay* , A^,A*(,, ...,Ap^, the 

minimization 

argmin f-/3^ - A^]3*^ + ^|^ - j3^|2 
Yk&kit V 2 

\ (34) 

+ E 

je.A'kUk ' ^ V 


is always attainable. 

Proof: First, note that the LMI conditions (l26l l- (l28T l are 
non-strict inequalities. The definition range of the solution 
matrices > 0,j € ^ are strict inequalities, but 

(l26l l- (l28T l imply that there exists a 5 > 0 such that X^>5 I 
and > 5 / for j G Thus, the feasible set Q,/; is closed 
and convex. 

Following again the proof of Theorem 1, (l26l l and (l28T l 
imply that are bounded for all k = 1,...,A^. In 

constrast, the variables Xj for j G are not restricted to a 
bounded set by the LMIs (I26l l- (l28l l. However, note that the 
cost function of (l34l l is quadratic in the variables Xj,j G 
Thus, due to the boundedness of we conclude that 

the sub-level sets of (|34]l 

{T,en,| -i5k_xkiik + ^\^_pk\2 

jeJiuk ^ ^ ^ } 

for c G M are bounded. Following the argument in Proposition 
4.1 in [9], Chapter 3, we can conclude that we can equiv¬ 
alently search for the minimum of the cost function over a 
non-empty sub-level set dTSl) instead of Hj.. Therefore, we 
can conclude that (l34l l is always attainable. 


Theorem 2: Algorithm 1 is a solution to Problem 1. 
In particular, the iteration steps Yi^{t),k = can be 

calculated in parallel, and satisfy the convergence conditions 

([I3]l, (0. 

Proof: Using Lemma [T] and |2] we can follow the steps 
from [9], Section 3.3 and 3.4, in order to prove convergence 
of the iterations. 


IV. Numerical Example 
Like in [6], we consider a system of the form ([T]), with 


A = 


B 


0.3775 0 0 

0.2959 0.3510 0 

1.4751 0.6232 1.0078 

0.2340 0 0 

0 0 0 

0 0 0 

[0.14 0] 


0 

0 

0 

0.5596 


0.4437 1.1878 0.0215 
0 2.2023 1.0039 

Di = 0.01 h for all k — 


which is observed by six sensor nodes, sensing two coor¬ 
dinates each. For every sensor an estimator is implemented, 
where none of the estimators is able to estimate the complete 
state vector without communication. The communication 
topology is assumed to be a directed circulant graph and we 
use Algorithm 1 to calculate the filter gains. For the numeri¬ 
cal calculations we use YALMIP [20]. Since we are dealing 
with a directed but balanced graph, we apply the method 
described in Remark 4 and use complete local representations 
of all variables Xj at every estimator k. The algorithm is 
run with both fixed performance parameter = 100 as 

discussed in Remark 3, and also using optimization over 
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Fig. 1. Evolution of the enor during iteration for fixed fJ. 
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Eig. 2. Evolution of the error (blue dots) and performance index (red 
line) in the algorithm involving optimisation over variables fJ*. 


the variable performance parameter (l29l l. 

where j3^(0) = 100 is set as initialization. 

In the first case, where j3 is fixed, we evaluate the 
matrix convergence condition (O by calculating the aver¬ 
age value Xj'’‘^ = subsequently Error = 

In the second case, involving optimisation over we 
additionally calculate ^ llk=i ^^d subsequently we 

have Error = + |j3^- 

The plots of the error evolution are shown in Figure [T] and 
|2] Figure |2] additionally shows the evolution of The 
graph demonstrates that is monotonically increasing, 
and since it is bounded from above according to Theorem 1, 
it must eventually converge to a limit. In fact, it eventually 
converges to 2.3 • 10^. 

Better performance j3 however is achieved at the expense 
of higher filter gains. For instance, after 70 iterations, the 
consensus gain Ki is 


21.1005 -0.0256 0.0196 
-0.0215 73.3369 0.5599 
-0.0423 -0.8806 99.8791 
-0.6033 -0.0178 0.0618 
0.0415 -0.0054 0.0972 
0.0117 -0.0003 0.0554 


-0.6018 0.0418 0.0117 
0.0073 0.0025 0.0021 
0.0617 0.0981 0.0536 
70.6692 1.3701 2.8005 
1.7726 20.7775 5.0466 
2.7740 3.2740 17.7281 








in the fixed-j3 case and 


28.8328 -0.0572 

0.0291 

-0.0089 

0.0333 

0.0655 

-0.0397 99.9887 

0.8962 

0.0198 

0.0050 

0.0044 

-0.0222 -0.8149 

100.0003 

0.2984 

0.1414 

0.0747 

-1.7921 -0.0044 

0.3121 

71.0412 

1.4502 

2.5260 

0.0751 -0.0035 

0.1384 

1.8242 

27.2461 

6.2623 

-0.0177 0.0016 

0.0711 

2.4613 

4.4543 

23.0606 


in the variable-j3 case. 


V. Conclusion 

We have developed a method for distributed filter design 
for cooperative .^^-type estimation. In order to achieve 
this we separated the centralized problem by introducing 
additional variables and then applied an algorithm that works 
locally and only needs communication for average consen¬ 
sus. 
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